Fifty years of data show the effects of climate on overall skull size and the extent of seasonal reversible skull size changes (Dehnel's phenomenon) in the common shrew

Abstract Global climate change affects many aspects of biology and has been shown to cause body size changes in animals. However, suitable datasets allowing the analysis of long‐term relationships between body size, climate, and its effects are rare. The size of the skull is often used as a proxy for overall body size. Skull size does not change much in fully grown vertebrates; however, some high‐metabolic small mammals shrink in winter and regrow in spring, including their skull and brain. This is thought to be a winter adaptation, as a smaller brain size reduces energy requirements. Climate could thus affect not only the overall size but also the pattern of the size change, that is, Dehnel's phenomenon, in these animals. We assessed the impact of the changes in climate on the overall skull size and the different stages of Dehnel's phenomenon in skulls of the common shrew, Sorex araneus, collected over 50 years in the Białowieża Forest, E Poland. Overall skull size decreased, along with increasing temperatures and decreasing soil moisture, which affected the availability of the shrews' main food source, earthworms. The skulls of males were larger than those of females, but the degree of the decrease in size did not differ between sexes. The magnitude of Dehnel's phenomenon increased over time, indicating an increasing selection pressure on animals in winter. Overall, climate clearly affected the common shrew's overall size as well as its seasonal size changes. With the current acceleration in climate change, the effects on the populations of this cold‐adapted species may be quite severe in a large part of its distribution range.


| INTRODUC TI ON
Global climate change affects animals in many ways, from phenology, geographic distribution, and phenotypic traits to distributions and population dynamics (Gardner et al., 2011;Parmesan, 2006;Walther et al., 2002). Changes in phenotypic traits correlated with global warming mainly involve temporal trends in body mass and size (Gardner et al., 2011), two traits that are pivotal for individual life histories (Sauer & Slade, 1987). Importantly, although mass and size are often used interchangeably, their change arises from different underlying physiological mechanisms. Mass is more sensitive to short-term environmental conditions affecting fat stores and the mass of the gastrointestinal tract and other organs (Canale et al., 2016;Hume et al., 2002;Piersma et al., 1999).

How climate influences body size is often described based on
Bergmann's rule, which in its commonly used version describes a positive relationship between body size and latitude of occurrence within species. Individuals of many homeothermic species are smaller at lower latitudes where temperatures are higher (e.g., Ashton, 2002;Ashton et al., 2000). This implies that a similar relationship should occur in response to temperature change over time. It has been postulated that a decrease in body size is the third universal response to contemporary global warming in addition to changes in distribution and phenology (Gardner et al., 2011). Although the majority of studies did not find significant temporal changes in the sizes of birds and mammals, there are many examples of long-term body size and body mass changes (recorded over one to several decades) correlated with global climate change. While a decrease in size is indeed the major response in birds, the situation is much less clear in mammals, where many species increase in size (Naya et al., 2017;Nengovhela et al., 2020;Teplitsky & Millien, 2014). This increase in size is usually explained as a reaction to increased food availability as a result of increasing temperatures (Boutin & Lane, 2014 and references cited therein;Yom-Tov & Geffen, 2011).
Soricine (red-toothed) shrews, especially Sorex shrews, are an excellent study system to investigate climate change-induced patterns in species-level body size changes. They have extremely high metabolic rates, much higher than expected for their body mass (Ochocińska & Taylor, 2005;Taylor, 1998). Consequently, shrews require a constant high food supply (Hanski, 1994;Keicher et al., 2017).
In addition, the body size of several Sorex species, represented by skull length, shows positive correlations with temperature and negative correlations with latitude across their distribution range, contrary to Bergmann's rule (Ochocińska & Taylor, 2003). It was suggested that food scarcity during winter in cold northern climates is a major factor selecting for small body size in shrews. This is in line with the "Resource Rule" of McNab (2010), which posits that "mammalian species will become larger or smaller depending on the size, abundance and availability of resources." A decrease in body size with latitude and in cold areas was also found in the masked shrew, Sorex cinereus, in Alaska and interpreted as related to better food resources in warmer areas (Yom-Tov & Yom-Tov, 2005). In agreement with this, the body size of S. cinereus in Alaska increased during the second half of the twentieth century along with increasing temperature, presumably due to increasing food availability in winter (Yom-Tov & Yom-Tov, 2005). However, an increasing temperature, especially when coupled with a decrease in precipitation, can also lead to a decrease in food resources. Drought has a negative impact on the abundance of shrews invertebrate prey (Coyle et al., 2017).
The unfavorable seasonal food conditions (caused directly or indirectly by climate) may bring about the advantage of being smaller in seasons when food is scarce. This is visible in the remarkable seasonal changes in size and morphology exhibited by soricine shrews.
The size of their skull and brain undergo a profound seasonal and reversible transformation: seasonal shrinkage and regrowth of individuals, whose magnitude is only comparable with that in some mustelid species LaPoint et al., 2017;Lazáro et al., 2017). Skull (braincase) height, but not skull length, decreases in winter, resulting in a flattened skull shape. This change in the skull, brain, and many other organs and tissues, as well as the body mass of soricine shrews and some other small high-metabolic animals, is called Dehnel's phenomenon Dehnel, 1949;Pucek, 1963Pucek, , 1970. After juveniles are fully grown, their brain mass decreases by 21% on average, and skull height decreases by 13% in S. araneus. Their brain mass reaches a minimum value in winter and then partially regrows in spring (Lazáro et al., 2021). Body mass, in contrast, decreases in anticipation of winter but then almost doubles in spring. Spring growth is related to reaching sexual maturity and likely to improved weather and food conditions. Sorex araneus, born in early summer, has a maximum lifespan of approximately 14 months, and almost all individuals die before the second winter after a terminal reproductive period (Pucek, 1981). The size decrease in winter is thought to be an adaptation to harsh climatic conditions in these nonhibernating, high-metabolic animals. The low body mass in shrews in winter has been hypothesized to reduce absolute food requirements when food availability is limited (McNab, 1991;Mezhzherin, 1964;Taylor et al., 2013). Accordingly, small winter animals have the same mass-corrected energy consumption as larger, first-and secondsummer shrews, even under ambient conditions with temperatures differing by as much as 30°C (Schaeffer et al., 2020). This results in absolute energy savings in winter and correlates with reduced food requirements. Reducing the size of the brain, an energetically expensive tissue (Aiello & Wheeler, 1995;Isler & van Schaik, 2006), may lead to further energy savings. Additionally, the decline in the relative weight of the liver, spleen, and adrenals and decreased synthesis by the endocrine glands in winter may contribute to lower metabolic rates in S. araneus and compensate for the increase in relative heart mass during this time (Hyvärinen, 1984;Pucek, 1965).
The magnitude of the winter decline in skull height in S. araneus, a characteristic component of Dehnel's phenomenon, increases toward the northeast in Europe and is linked to large-scale environmental conditions and probably also the local habitat structure (Lazáro et al., 2021;Pucek, 1970). It is positively correlated with temperature seasonality, annual temperature range, and other climate parameters, although no such relationships have been observed with skull height regrowth (Lazáro et al., 2021), which has led to the hypothesis that shrinking may be the result of different evolutionary drivers than regrowth. In line with these observations, the individual decrease in skull height in S. araneus is flexible and modulated by ambient temperature (Lazáro et al., 2019).
Our aim was to reject a null hypothesis that the change in climate over several decades has had no impact on (1) the overall size of S. araneus as measured by skull dimensions and (2) the intensity of the reversible seasonal size change (Dehnel's phenomenon We first hypothesized that H1: because they are larger at lower, warmer latitudes, the body size of S. araneus increased over the years as a response to increasing temperatures (and increasing food availability as in the case of S. cinereus) or Ha1: decreased as a response to a drier climate leading to a decreasing soil water level and thus lower earthworm availability in the Białowieża Forest. We further hypothesized that H2: increasing winter temperatures reduce the winter decrease in skull height (and thus brain size). Conversely, Ha2: the winter decrease in skull height will become greater if climate change is combined with increasingly less favorable food availability in winter than in summer. Academy of Sciences, in Białowieża. Shrews were trapped in a mesic broadleaved forest with a transition to a moist broadleaved forest, mainly in compartment No. 371 (52°44′N, 23°52′E) but also in adjacent Nos. 369 and 370 (see Dehnel, 1949 for a map; all these compartments are in the strict reserve). This type of forest is classified as a subcontinental oak-lime-hornbeam forest Tilio Carpinetum. The mean age of the trees exceeded 100 years, and the amount of dead wood on the ground averaged 50 m 3 ha −1 .

| Shrews, study area and trapping
Sorex araneus were caught with pitfall traps (cylinders 30 cm deep and cones 45 cm deep). After 1965, wooden box live traps were additionally used. We selected and measured the skulls of shrews from the three main stages of Dehnel's phenomenon: sexually immature summer individuals in their first calendar year (hereafter "juveniles") caught from July to August; smaller, immature, winter individuals ("subadults") caught from December to February; and regrown, sexually mature spring/summer individuals in their second calendar year ("adults") caught from April to August. The individuals were assigned to the above three age classes upon capture, and their age was recorded in the collection database. Their sex was determined by dissection. Juveniles are easily distinguishable from adults by the degree of tooth wear, development of gonads, and fur appearance (Churchfield, 1990;Pankakoski, 1989). Because S. araneus lives for approximately 1 year and only reproduces in spring and summer, there is no generation overlap of mature adults. In winter, there is only one age class-subadults. After 1981, juvenile individuals were trapped only in July. There was no winter trapping at all after 1985.
We also measured three individuals each from June (juveniles) and November (subadults) and two from September (juveniles) from 1955 to 1985.

| Measurements of skull dimensions
A single observer measured the skull height, length, and width (mm) using X-ray images of the alcohol-preserved specimens, as described by Lazáro et al. (2017). Briefly, skull length was measured from the anterior-most projection of the first incisive tooth to the occipital condyle ( Figure  the % of the mean) for the skull length, height, and width was 1.35%, 2.04%, and 1.89%, respectively (supplemental information in Lazáro et al., 2017). We measured a total of 502 skulls.

| The choice of months representing the three age categories
We compared skulls in three age categories of individuals who characterize their entire lifetime: fully grown juveniles, size-decreased F I G U R E 1 The common shrew (Sorex araneus). Photo by Leszek Rychlik.
subadults, and regrown adults ( Figure 3). To test the homogeneity of each age category, we used an ANCOVA model with the month as a factor, the time (year) as a covariate, and the interaction between month and year. The skull height of juvenile individuals differed between July and August (F 1,201 = 8.08, p = .0049), and the interaction was significant (F 1,201 = 8.12, p = .0048) in the period between 1953 and 1981, when individuals from both months were measured. Thus, we further analyzed July and August juveniles separately.
We chose the measurements of the skull heights from January to February, when they are the lowest during the lifespan of S. araneus We ran the same tests for skull width and length, which vary much less between seasons (Lazáro et al., 2017). The year × month interaction was not significant in any age category, and there were no significant differences between collection months. Including sex in the above models did not change the results. Table 1 1978, 1979, 1983, and 1987, when shrews were caught in the second half of May, but not in June, no juvenile S. araneus were caught, only adults were caught. In 1981, when no catching was performed in May, the first juvenile was caught on June 1. In summary, we found no signs of change in the breeding phenology, at least between 1953 and 1990. Moreover, reaching the minimum skull height in winter did not change over the years; it occurred in January and February, regardless of the decade or locality (Dehnel, 1949;Lazaro et al., 2019;Pucek, 1970).

| The change in skull size from 1953 to 2004
We first analyzed the changes in the overall size of S. araneus using the full dataset from 1953 to 2004, which was only available for juveniles from July and adults. We tested the temporal change in all three skull dimensions, skull height, width, and length, with three F I G U R E 2 X-ray images of Sorex araneus skulls and linear measurements taken from the images: (a) skull length (SKL), skull (braincase) height (BCH), and (b) skull width (BCW). See text for details.

F I G U R E 3
Skull height throughout the lifetime of Sorex araneus from Białowieża Forest. Sampling months and corresponding age categories are shown. The first weaned shrews appeared in June (juveniles); they then shrank toward winter (subadults) and finally regrew in the next spring/summer (adults), when the adults reproduced. Most of the shrews then died before the next winter. Solid circles are the means ± standard errors when the sample size was ≥5, open circles are the means with sample size <5. Horizontal bars represent age categories and months from which skulls were analyzed; skulls of juvenile shrews from July and August were analyzed separately (see text). The presented data were limited to the period between 1955 and 1985, for which we had individuals from all three age categories.
ANCOVA models with age and sex as factors and time (year) as a covariate. We added a quadratic term to the linear models to test for linearity of the relationships and determine whether it improved the model according to an F test with a type I sum of squares.
We also tested the linearity of the regressions of the three dimensions over time for two age categories separately. The quadratic term was only significant for the skull length of adults (p = .032).
However, inspection of the residuals from other linear regressions suggested a more subtle nonlinearity of the skull size over time.
Thus, we additionally used locally weighed regression smoothing (LOESS; Cleveland, 1979). We selected the smoothing parameter by minimizing Akaike's information criterion, which strikes a balance between the residual sum of squares and the complexity of the fit.

| The influence of climate change on the size of S. araneus
We examined the weather data from Białowieża from 1952 to 2004 to explain changes in skull size in juvenile and adult S. araneus. We obtained the mean daily temperature (°C), rainfall (mm), and snow cover (mm) from the Białowieża Meteorological Station. We summed or averaged these daily measures for longer periods and computed the soil moisture deficit from monthly temperatures and precipitation using the Watbug program (Willmott, 1977). Monthly deficits were summed for years, or other periods as indicated below. We used these deficits as the indicators of earthworm abundance in the soil. As in the previous section, we used linear regression, tested for nonlinearity by adding quadratic terms, and smoothed the relationships with LOESS.

F I G U R E 4
The relationships between (a) skull height, (b) skull length and year in July juveniles (JUV-7;1953-2004, August juveniles (JUV-8;1953-19811955-1985 and adults (AD;1953-2004 of Sorex araneus compared with changing climate parameters: (c) the mean annual temperature, (d) the annual sum of precipitation, and (e) the annual sum of soil moisture deficit. The 95% confidence limits of the relationships smoothed with LOESS are shown (for the sake of clarity, confidence limits for JUV-8 and SAD in (b) are not presented). Additionally, for clarity, both sexes are combined in (a) and (b) (the difference in skull size between sexes was small and the change with time was parallel; see text).
We calculated the mean daily temperature, sum of the precipitation, and sum of the moisture deficit in the soil in the period preceding the collection of specimens (from their approximate time of birth; see Section 3 for details). We also used the mean temperature of January and February, the number of days with snow cover in these months and during the whole winter preceding the collection of adult S. araneus.
To select models best characterizing the relationship between climatic variables and skull dimensions we ordered all possible models (separately for juveniles and adults) according to the AICc (sample size corrected AIC). Then we inspected all models with delta <3 for most influential variables. Next instead of averaging coefficients for equally supported models (as for all skull measurements within the two age classes there was no single best model), we selected climatic variables which are not correlated and significantly explain most of variation (based on R 2 ). We checked for a temporal autocorrelation (correlation of residuals in neighboring years) for all three skull measurements using autocorrelation function (ACF) in R for a wide range of lags. Autocorrelation was not significant in any of the models.

| Skull dimensions in all three age categories from 1955 to 1985
We compared absolute values of the skull height, width, and length for the years when specimens for all three age categories were available in an ANCOVA with age and sex as the factors and sampling year as the covariate.

| Change in the magnitude of Dehnel's phenomenon from 1955 to 1985
We compared the amount of change in skull height from the subset where we had data from all three size extremes (summer juveniles, winter subadults, and regrown adults). To do this, we compared juveniles (July and August separately) with subadults and subadults with adults in ANCOVA models with age as the factor, year as the covariate, and the interaction between year and age.
To visualize the change in the magnitude of Dehnel's phenomenon over time, we divided the years where we had data from all three size extremes (summer juveniles, winter subadults, and regrown adults) into three time intervals (1955-1964, 1965-1974, and 1975-1985). We divided skull height by skull length (which does not change within individuals) to control for differences in overall body size between individuals. Next, we compared the seasonal curves of the relative skull height between these three time intervals with a generalized additive model (GAM) using "time interval" as a parametric term and "month" as a nonparametric term. We used a Gaussian distribution and applied a smoothing function to "month" with knots restricted to five (because we had an a priori expectation for the presence of the shape of a seasonal skull height pattern).

| Statistics
We used general linear models (GLMs) for ANOVA and ANCOVA  Figure 4a,b; see Figure S1 for scatterplots and regression lines). Skull height and length, but not skull width, differed between the two age categories and were larger in juveniles than in adults (Table 2a,b).
The interactions between sex and time (years) were not significant when included in the models (F 1,343 = 0.55, p = .458, skull height; F 1,343 = 0.37, p = .543, width; F 1,340 = 0.28, p = .598, length), which indicated that male and female skulls were shrinking at the same rate. All three measures were larger in males than in females, except for adult skull length (Table 2a,b, Figure S1). The largest difference between the two sexes was in skull height: 2.8% in adults and 1.1% in juveniles (Table 2b). Characteristically, all relationships of skull dimensions with time (except skull height in adults) revealed a common pattern when smoothed with LOESS: little change until the early 1970s and decreases thereafter (Figure 4 and Figure S2).

| The influence of climate change on the size of S. araneus
The weather records from this period revealed several significant trends indicative of a changing climate (scatterplots of six climate parameters are shown in Figure S3). The mean annual temperature increased significantly (F 1,51 = 28.1, p < .0001; Figure 4c and Figure S3a) from 6.1 to 8.2°C (as calculated from the regression line).
There was also an increase in the mean July temperature (F 1,51 = 10.5, p = .0021; Figure S3b We found that soil moisture deficit and precipitation were correlated with skull size in S. araneus and could be the drivers of skull size decrease over years ( Table 3). Although skull length, the proxy of the overall body size, was not correlated with climate parame-  TA B L E 2 (a) Results of three ANCOVA models testing the relationship between skull height, width, and length and age category (juveniles from July and adults), sex, and year of sampling in Sorex araneus in Białowieża Forest from 1953 to 2004; n = 348 individuals (n = 345 in skull length; without three outliers in juveniles). All significant effects are marked in bold. All four possible interactions added to each of three models were nonsignificant (b) Mean dimensions of skulls (mm) with standard errors for juvenile and adult males and females; different letters (a, b) mark significant differences of a given dimension between the sexes in the same age category.

TA B L E 3
The relationship of the change in skull height, length, and width in juvenile (July) and adult Sorex araneus with weather variables between 1953 and 2004, and in skull length between 1971 and 2004 (the period of significant increase in drought). Note: SumDef, T, and SumPrec are shrew's "whole-life" weather parameters; for juveniles, they were calculated for the period between May and June, which covers the early life of juveniles and preceded their collection in July; for adults, they were calculated from June, when the first cohort of shrews was weaned, to March of the following year, when adult shrews were collected between April and August. SnowJF = the number of days with a snow cover in January and February preceding the collection of adults. R 2 is the proportion of the variance for the skull size explained by the independent variables in the multiple regression. Significance of the independent variables: * = 0.05; ** = 0.01; *** = 0.001.

Model
Abbreviations: SumDef, the sum of monthly soil moisture deficit; SumPrec, the sum of precipitation; T, the mean temperature.
found an impact of winter weather (during the subadult stage) on adult skull height the following summer: summer skull heights were greater when there were more days with snow cover in the previous January and February (Table 3).

| Skull dimensions in all three age categories from 1955 to 1985
In the time period from 1955 to 1985, we had measurements from all three ages: juveniles, subadults, and adults. Skull height, width, and length were significantly correlated with age (confirming the typical seasonal size change pattern; Table 4a). Only skull height showed a negative significant relationship with year, and the interaction between year and age was significant when included in this model.
All measures of skull size differed significantly between the age categories (Table 4b). Skull height decreased significantly by 14.4% from summer juveniles to winter subadults and regrew by 7.4% in the following summer adults. Skull width followed this pattern with a lesser but statistically significant decrease of 1.4% and then an increase of 1.9%. The 2% decrease in skull length from juveniles to subadults was also statistically significant (Table 4b).

| Change in the magnitude of Dehnel's phenomenon from 1955 to 1985
We compared the slopes of the regression lines of skull height over time (years) between juveniles versus subadults and subadults versus adults (Figure 4a, see Figure S4 for regression lines and scatterplots).
The skull height of July juvenile shrews did not change over the tested period (F 1,149 = 1.34, p = .249) but decreased significantly in winter subadults (F 1,68 = 19.52, p < .0001; Figure 4a), which was reflected in a significant interaction between age and time (years) (F 1,217 = 8.23, p = .0045). This means that the decline in skull height from July to winter was greater with time. The interaction was not significant when F I G U R E 5 Relationship between the skull length in adult Sorex araneus and the soil moisture deficit; regression line with the 95% confidence limits; the lowest, outlying point of skull length was omitted from the regression. The period for which the soil moisture deficit was calculated is explained in Table 3.  the regression lines of August juveniles and subadults were compared (F 1,118 = 0.63, p = .430). The skull height of adults did not change with time (F 1,72 = 0.00, p = .985); considering the steep decrease in subadults, the regrowth in skull height from subadults to adults increased over time (interaction: F 1,140 = 9.54, p = .0024; Figure 4a). Sex, when added to these models was not significant.
This phenomenon is also illustrated in Figure 6, which depicts the change in shrew skull height corrected for skull length during its lifetime over three decades. The standard errors of juvenile and adult skull height overlapped considerably, while subadult skull height decreased strongly over time (years). In summary, the magnitude of Dehnel's phenomenon increased over the tested period due to winter animals becoming smaller.

| Decrease in skull and thus overall size over time (1953-2004)
Skull dimensions, skull length, the often-used proxy for overall body size, as well as skull height and width, decreased significantly in juvenile and adult S. araneus during the study period. We found sexual dimorphism in skull size; males had larger skulls than females, but the sexes did not differ in the degree of decrease in skull size. Our results add to the very mixed picture of the change in body size in mammals under global climate change. The analysis of 50 species of rodents from seven families revealed no significant temporal change in body size (represented mainly by skull dimensions) in 29 species, a significant decline in 13 (26% of the total number), and an increase in size in eight species (16% of the total) (Nengovhela et al., 2020). Of the three carnivore species in which skull size was measured in the second half of the 20th century and analyzed in the context of climate change, skull size showed a significant increase in two species and a significant decline in the other (Yom-Tov et al., 2008;Yom-Tov, Leader, et al., 2010;Yom-Tov, Roos, et al., 2010). Thus, S. araneus in Białowieża adds to the number of rather rare cases of decreased body sizes among mammals. However, when body mass is analyzed in terrestrial mammals, a clear tendency of the decrease in body mass is observed; this is also true for shrews in the 20th century (Naya et al., 2017).
The decrease in body size in S. araneus in Białowieża is in contrast with the observed increased body size (body and foot length) in S. cinereus in Alaska during the second half of the twentieth century along with increasing temperature. In S. cinereus, the increase was presumably due to higher food availability in winter as a result of higher winter temperatures, favorable for its prey (Yom-Tov & Yom-Tov, 2005), which consists mainly of insect larvae, beetles and spiders (Whitaker, 2004). We argue that the smaller size of S. araneus in summer could be an effect of increasing drought leading to reduced availability of earthworms, its main food source. The negative impact of decreased soil moisture on the abundance of epigeic and endogeic invertebrates, especially earthworms, has been well documented (Coyle et al., 2017;Singh et al., 2019). Earthworms have very limited morphological and physiological means for reducing water loss through the cuticle and need moist soil to avoid dehydration (Carley, 1978). Sorex araneus strongly depends on a high abundance of food because its metabolic rates measured under natural conditions equal 258% of the value predicted for an average mammal of the same body mass (Ochocińska & Taylor, 2005). Consequently, S. araneus prefers humid habitats with humus-rich soil (Shchipanov et al., 2019), such as our study plot, which harbor large numbers of earthworms.
In support of this, the decrease in skull length in adults be-  Table 3 also support the importance of these two climate parameters for the long-term decrease in body size of

S. araneus.
There were no regular measurements of the soil water level in the Białowieża Forest before 1985, and those that existed were not comparable with later measurements. However, between 1985 and 2001, the water level decreased in mixed deciduous biotopes by 40 cm. This resulted from an increase in temperature, a decrease in rainfall, and winters without snow and thus a lack of water accumulation at the beginning of the vegetation period (Pierzgalski et al., 2002). The progressive decrease in the soil water level was certainly magnified by the intense drainage work on peatlands in the eastern Belorussian part of the Białowieża Forest in the late 1950s.
The average biomass of earthworms in the upper soil layer in the mixed deciduous areas of the Białowieża Forest was 42.2 g m −2 between 1997 and 2000 and depended on soil moisture: it was highest in spring, decreased to 22% of the spring values in summer, and increased to 41% in autumn (Kowalczyk et al., 2003). This pattern followed the annual course of soil humidity in the Białowieża mixed deciduous forest (Pierzgalski et al., 2002). Earthworms were not available or were hardly available in the forest in winter (Kowalczyk et al., 2003).  (Kowalczyk et al., 2003). The shift to prey other than earthworms, the reduction in the population density of S. araneus, and the observed decrease in body size (and amount of consumed food) may explain why no change in the breeding phenology of S. araneus was detected (see Section 2) despite climate change.
While the increasing soil moisture deficit and decreasing precipitation affected both juveniles and adults, resulting in smaller skull size, the adults were additionally affected by the conditions in the preceding winter. Milder minimum temperatures in February and a decreasing number of days with snow cover in the coldest months of January and February resulted in adults with lower skull heights ( Table 3). The lack of snow cover leads to freezing of the upper soil layer and substantially reduces the abundance of various arthropods, which are important components of the winter diet in S. araneus (Templer et al., 2012). Snow generates a relatively mild and stable subnivean microclimate in the litter layer where shrews live (Churchfield et al., 2012;Coulianos & Johnels, 1962), which is not only favorable for their prey but may also diminish the costs of thermoregulation and consequently the food requirements of the shrews (Aitchison, 1987).  (Ochocińska & Taylor, 2005), fits this pattern. It is likely that the relationship between body size decrease and high metabolic rate is generally true for soricine shrews. Three out of four American Sorex species (not including S. cinereus) had a decreased body mass during the 20th century, and two of them had a significantly decreased body mass (Naya et al., 2017). The body mass of the short-tailed shrew (Blarina brevicauda) also decreased in 20th century (Naya et al., 2017). However, data from different ages and thus Dehnel stages were pooled in this study, potentially affecting the results and hiding the effects of seasonal body size changes.

| Changes in the magnitude of Dehnel's phenomenon (1955-1985)
We found that the decline in skull height from summer to winter in S. araneus from the Białowieża Forest became stronger between 1955 and 1985 ( Figure 6). We standardized skull height with skull length in the analysis of Dehnel's phenomenon to describe changes in relative skull height. Relative skull height is strongly correlated with brain size (Bielak & Pucek, 1960;Lazáro et al., 2018). Thus, our results also reflected changes in the magnitude of brain size changes.
Initially, the skull height of juveniles did not differ between July and August from 1955 to the 1970s. The skull height in these 2 months then diverged over the following years ( Figure 4a). The smaller skull height in August than in July after 1970 meant that the decrease in skull height started earlier, already in summer. The magnitude of the decrease was also higher because the relative skull height of winter subadults decreased even more over time ( Figure 6). We also found a significantly greater skull height regrowth from 1955 to 1985, which resulted from a steep decrease in skull height in the subadults, while adults did not change (Figure 4a). Almost nothing is known about the drivers of this regrowth, which are presumably different from those of the decrease and are likely associated with preparation for the terminal reproductive period in shrews (Lazáro et al., 2019(Lazáro et al., , 2021. Comparison of the skull height regrowth in different S. araneus populations revealed only one correlation with weather parameters, namely, a positive correlation with precipitation during the warmest quarter of the year (Lazáro et al., 2021). This also points to changes in food availability.
In summary, we documented a decrease in overall size as represented by absolute skull height and length in the 52-year  skull series and identified increasing temperature and drought as the main correlates of the decrease. Interpreting the results regarding Dehnel's phenomenon, which increased in magnitude because of a sharp size decrease in winter subadults, was more difficult. Our dataset about winter subadults stopped shortly before the time point when climate change became significant, and it would have been interesting to see how subadult size developed after this.
There are strong indications that there are important effects of climate on food availability that indirectly cause trends in both skull size and the magnitude of Dehnel's phenomenon. However, this must be clarified experimentally.
With continued climate change, all the parameters that we found to be of importance for the size of the common shrew and its unique winter adaptation are expected to increase even more strongly. This may negatively affect populations of S. araneus in many regions that may become too hot and/or dry for these cold-adapted animals.  commented on the previous version of the manuscript. We thank two anonymous reviewers for their constructive comments. We did not receive any grants to conduct this research.

CO N FLI C T O F I NTE R E S T
The authors declare that they have no conflicts of interest.

DATA AVA I L A B I L I T Y S TAT E M E N T
The primary data on skull sizes and climate in Białowieża can be found in Dryad: https://doi.org/10.5061/dryad.xd254 7dm7.